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I.  INTRODUCTION 


There  are  many  classical,  nonlinear  systems  exhibiting  some  kind  of 
chaotic  behavior.  Examples  include  the  turbulent  flow  of  a  fluid, 
usually  described  by  means  of  the  Navier-Stokes  equations,  and  the  be¬ 
havior  of  liquids,  gases  or  antiferromagnets  above  the  critical  point, 
among  others. 

In  this  paper,  we  propose  to  reexamine  and  further  develop  an  ap¬ 
proach  to  the  description  of  such  systems,  originally  proposed  by 

1  2 
Martin,  Siggia  and  Rose,  and  reformulated  by  DeDominicis  and  Peliti. 

This  approach  promises  to  yield  some  qualitative  insight  into  problems 

of  the  kind  outlined  above;  moreover,  if  sufficiently  developed,  it  may 

be  combined  with  numerical  methods  in  order  to  increase  the  accuracy  of 

the  results  or  reduce  computing  time.  We  use  the  formulation  of  ref.  2 

and  restrict  our  attention  to  autonomous  systems.  From  a  qualitative 

point  of  view,  the  method  consists  of  the  following.  One  considers  a 

system  described  by  a  set  of  deterministic  nonlinear  equations  governing 

the  behavior  of  the  dynamical  variables.  In  the  approach  of  refs.  1  and 

2,  the  system  is  then  perturbed  by  a  random  force.  It  is  then  shown 

that  this  theory  of  a  "randomly  stirred"  system  is  formally  equivalent 

to  a  quantum  field  theory:  one  can  define  it,  for  instance,  in  terms  of 

a  suitable  path  integral,  and  compute  fluctuations  induced  by  the  random 

stirring  force  by  means  of  familiar  and  internally  consistent  methods. 

However,  the  approach  as  described  in  refs.  1  and  2,  is,  to  a 
certain  degree,  incomplete  from  a  physical  point  of  view.  Jt  is  not 
clear  that  a  system  has  to  be  randomly  stirred  in  order  to  become 
chaotic:  for  example,  a  fluid  flow  can  apparently  become  turbulent 

under  the  influence  of  some  perfectly  smooth,  non-fluctuating  external 
force,  such  as  a  constant  pressure  gradient  in  a  pipe.  Although  one  may 
argue  that  in  a  "real"  flow  there  is  always  some  randomness  present 
(such  as  thermal  fluctuations  in  the  fluid,  roughness  of  the  wall  of  the 

pipe,  etc.),  a  recent  numerical  simulation  carried  out  by  Deissler3 
seems  to  indicate  that  this  is  not  crucial  in  the  development  of  turbu¬ 
lence  and  probably  in  the  chaotic  behavior  of  other  systems  either. 
Intuitively,  one  would  like  to  "switch  off"  the  random  stirring  force  in 
some  smooth  manner,  and  make  sure  that  chaotic  behavior  persists  even  if 
the  random  perturbation  is  weak.  (This,  to  a  certain  extent,  would 


*P.  C.  Martin,  E.'D.  Siggia  and  H.  A.  Rose,  "Statistical  Dynamics  of 
Classical  Systems,"  Phvs ,  Rev. ,  Vol.  A8,  pp.  423-440,  1973. 

2 

C.  DeDominicis  and  L.  Peliti,  "Field-Theory,  Renormalization  and  Criti¬ 
cal  Dynamics  Above  T  :  Helium,  Antiferromagnets  and  Liquid-Gas 

c 

Systems,"  Phvs.  Rev. ,  Vol.  B18,  pp.  353-376,  1978. 

3 

R.  G.  Deissler,  "Is  Navier-Stokes  Turbulence  Chaotic?,"  Phvs.  Fluids.. 
Vol.  29,  pp.  1453-1457,  1986. 
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justify  the  conjecture  expressed  by  Martin  and  DeDorainicis  that  at 
least  some  important  properties  of  a  turbulent  flow  are  independent  of 
the  agitation  mechanism.)  From  the  technical  point  of  view,  the  formal¬ 
ism,  as  it  stands,  does  not  readily  lend  itself  to  the  development  of  a 
systematic  loop  expansion  around  some  non-trivial,  average  solution 
represented  by  a  non-trivial  saddle  point  of  the  path  integral,  unless 
the  stirring  force  is  kept  nonzero,  since  all  propagators  are  propor¬ 
tional  to  the  correlation  function  of  the  stirring  force. 

Quite  clearly,  this  calls  for  a  more  careful  examination  of  the 
role  played  by  a  random,  typically  Gaussian,  stirring  force  in  the 
theory  of  nonlinear  systems  in  which  we  are  interested.  It  turns  out 
that,  within  the  framework  of  the  formalism  created  by  Martin,  Siggia, 
Rose,  DeDominicis,  and  Peliti  (MSRDP) ,  loc.  cit. .  one  can  associate  a 
Hamiltonian  system  with  every  autonomous,  but  possibly  dissipative 
system  stirred  by  a  random  force.  That  Hamiltonian  system  becomes 
singular,  in  a  sense  to  be  specified  later,  if  from  the  outset,  the 
stirring  force  is  set  equal  to  zero.  In  some  simple  and  physically 
interesting  cases,  typically,  Gaussian  stirring  of  equal  strength  at 
every  frequency  and  wave  number,  the  functional  integration  over 
canonical  momenta  can  be  carried  out  in  a  closed  form.  This  establishes 
a  relationship  between  the  Hamiltonian  and  Lagrangian  forms  of  the  path 
integral,  as  long  as  the  stirring  force  is  not  exactly  zero. 

The  paper  is  organized  as  follows.  In  the  next  Section  we  consider 
some  general  properties  of  systems  described  by  the  class  of  equations, 
of  the  form  c^X  -  F(X)  =  ef(x,t),  where  f  is  a  Gaussian  stirring  force 

and  e  is  a  parameter  characterizing  the  strength  of  its  coupling  to  the 
system  of  interest,  while  all  dynamical  variables  are  collectively  de¬ 
noted  by  X,  where  X  is  an  element  of  some  vector  field.  We  explicitly 
exhibit  the  correspondence  between  the  MSRDP  path  integral  and  a 
Hamiltonian  system  and  derive  the  Lagrangian  form  of  the  path  integral 
which  is  suitable  for  the  development  of  a  systematic  approximation 
scheme.  In  Section  3  the  formalism  is  applied  to  the  turbulent  flow  of 
a  viscous  incompressible  fluid.  In  particular,  we  derive  the  general 
form  of  the  Feynman  rules  necessary  for  the  computation  of  the  non- 
Gaussian  part  of  the  velocity  distribution.  In  Section  4  this  formalism 
is  illustrated  on  a  flow  which  is  Poiseuille  on  the  average;  in  partic¬ 
ular,  we  compute  the  two-point  velocity  correlation  function  in  the 
Gaussian  approximation.  The  results  are  summarized  and  discussed  in 
Section  5. 


2.  GENERAL  FORMALISM 

For  the  sake  of  completeness,  we  start  with  a  brief  review  of  the 
procedure  of  DeDominicis  and  Peliti  (ref.  2).  Consider  a  "classical" 
system  described  by  an  equation  of  the  type 


4 

P.  C.  Martin  and  C.  DeDominicis,  "The  Long  Distance  Behavior  of 
Randomly  Stirred  Fluids,"  Prog.  Theor.  Phvs.  Sudd. .  No.  64,  pp.  108- 
123,  1978. 
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*tX  -  F (X)  =  ef, 


(2.1) 


where  f  is  a  Gaussian  random  force  and  t  is  a  numerical  parameter,  which 
ultimately  will  be  made  vanishingly  small.  The  correlation  operator  of 
f  is  denoted  by  K  and  the  scalar  product  is  then  denoted  by  the  symbol 
<  ,  >• 

The  principal  result  of  DeDominicis  and  Peliti  is  that  the  charac¬ 
teristic  functional  of  the  stochastic  process  defined  in  the  preceding 
can  be  written  in  the  form: 


G  [  j  ]  =  Jdx  Df  Det(<5tH5F/5X)6(dtX-F(X)-fef  )e1<j,X>e  <f,K  f  >/2  .  (2.2) 

Here,  DX  and  Df  stand  for  an  appropriately  normalized  functional  measure 
and  j  is  an  arbitrary  source  function.  The  cumulants  of  the  distribu¬ 
tion  are  generated  by  the  functional  W[j]  =  i  1lnG[j].  In  what  follows, 
we  are  mainly  interested  in  the  cumulants  of  any  given  distribution. 
Hence,  it  is  convenient  to  define  an  equivalence  class  of  functionals  of 
the  source  j  by  calling  two  functionals  equivalent  if  their  ratio  is 
independent  of  the  source  j.  Following  the  tradition  of  quantum  field 
theory,  members  of  this  equivalence  class  are  identified  by  the  conven¬ 
tional  equality  sign  (=).  With  this,  we  now  have,  using  a  Fourier 
representation  of  the  delta  functional , 


G[j] 


DX  Df  DP  Det(dt-5F/5X)e 


i< j , X>  ~<f , K 


-1 


f>/2ei<P, (d  X-F-6 f ) > 


(2.3) 


Integration  over  the  stirring  force  can  be  readily  performed.  The 
result  is: 

G[  j  ]  =  fl)X  DP  e1<:i  ,X>Det(at-6F/5X)exp  i{<P,3tX>  -  H),  (2.4) 


2 

where  H  =  <P,F>  +  ie  <P,KP>/2.  (We  note  that  the  correlation  operator 
is,  by  necessity,  a  Hermitean  one,  thus,  <f,Kg>  =  <Kf,g>  for  any  admis¬ 
sible  pair  of  vectors  f  and  g. )  In  what  follows,  one  can  set  the  cor¬ 
relation  operator,  K,  proportional  to  the  unit  operator  without  any  sub¬ 
stantial  loss  of  generality.  (Indeed,  such  a  choice  of  K  is  desirable 
from  an  intuitive  point  of  view:  one  "stirs"  the  system  by  a  white 
noise,  i.e.,  the  correlation  function  of  the  stirring  force  is  propor¬ 
tional  to  6(t-t ' )  and  the  stirring  is  of  the  same  strength  at.  every  wave 
number.)  With  this  simplification,  let  us  now  take  a  look  at  the  func¬ 
tional  integral,  eq.  (2.4).  Quite  obviously,  this  integral  is  of  the 

canonical  form  (see  e.g.,  Itzykson  and  Zuber'  ).  One  is  integrating  over 
the  phase  space  of  variables  X,P,  with  a  weight  exp  i{<P,d  X>  -  H). 

Clearly,  the  quantity  in  curly  brackets  can  be  identified  with  the 


5 

C.  Itzykson  and  J.  B.  Zuber,  " Functional  Methods,"  Quantum  Field 
Theory .  McGraw-Hill,  New  York,  1980,  Chapter  IX,  pp.  425-474, 
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canonical  form  of  the  action  governing  the  fluctuations,  while  H  plays 
the  role  of  the  Hamiltonian.  The  important  point  to  notice  is  that  if 
one  insisted  on  discarding  the  coupling  between  the  stirring  and  the 
system  under  consideration,  the  latter  would  become  a  singular  Hamilton¬ 
ian  system  (in  the  sense  of  classical  mechanics):  one  cannot  pass  free¬ 
ly  between  the  Hamiltonian  and  Lagrangian  formulations  of  the  problem  at 

2 

hand  since  the  second  derivative  operator,  6  H/5P6P,  is  singular  (in 
this  case,  identically  zero)  when  e  =  0.  (Given  the  fact  that  the  inte¬ 
gration  over  the  canonical  momenta,  P,  is  a  Gaussian  one,  reflecting  the 
Gaussian  nature  of  the  stirring,  there  is  a  one-to-one  correspondence 
between  the  analysis  of  the  "classical"  Hamiltonian  system  defined 
implicitly  by  eqs .  (2.3)  and  (2.4)  and  the  actual  system  under  consider¬ 
ation  .  ) 


Let  us  now  accept  a  unit  correlation  operator  of  the  stirring 
force.  One  can  then  readily  perform  the  integration  over  the  canonical 
momenta,  thus  arriving  at  the  "Lagrangian"  form  of  the  functional  inte¬ 
gral,  viz.  , 


G[j]  - 


DXe1<^ ,X>Det(dt-6F/6X)exp-e  2l/2<(atX-F) , (3tX-F)> . 


(2.5) 


In  general,  it  is  possible  to  exponentiate  the  functional  determinant, 
appearing  in  eqs.  (2.4)  and  (2.5)  by  introducing  the  usual  Fadeev-Popov 
ghost  fields,  cf . ,  Itzykson  and  Zuber,  lot:,  cit..  We  denote  the  ghost 
fields  by  u  and  v;  then  using  the  standard  expression  of  a  functional 
determinant  in  terms  of  an  integral  over  Grassmann  fields,  we  get  the 
final  expression  for  the  generating  functional: 

G [ j ]  =  jDX  Du  Dv  e1<j,X>exp  -  (t ~2 l/2< (3 tX-F ) , (3 fcX-F) >  -T),  (2.6) 

where  f  stands  for  the  ghost  part  of  the  action,  viz., 

r  =  <u ,  (d  -  6F/6X) v> .  (2.7) 

The  expression  of  the  characteristic  functional  obtained  in  eq.  (2.6) 
and  (2.7)  has  the  advantage  that  it  allows  one  to  perform  a  consistent 

2 

loop  expansion  (in  powers  of  e  )  in  a  standard  fashion,  cf.  Itzykson  and 

2 

Zuber,  loc ,  cit . ,  Ch.  6,  In  essence,  e  plays  the  same  role  in  this 
formalism  as  Planck’s  constant  does  in  quantum  field  theory.  It  fol¬ 
lows  that  it  is  impossible  to  put  e  =  0  from  the  outset:  this  is  a 
singular  limit  in  the  same  sense  as  the  limit  h  -»  0  is  in  quantum 
theory.  In  a  general  system  all  one  can  do  is  to  rescale  the  variables 
of  the  functional  integration  in  a  manner  analogous  to  a  quantum  field 
theory  with  a  single  coupling  constant;  in  this  way,  a  "small  e,"  expan- 
sion  can  be  converted  to  a  "weak  nonlinearity"  expansion.  A  practically 
important  exception  to  this  scaling  argument  is  given  by  a  scale  invar¬ 
iant  system,  notably  the  hydrodynamics  of  a  viscous  incompressible 
fluid.  In  this  case,  a  scale  transformation  can  be  chosen  such  that  the 
small  parameter  e  is  altogether  eliminated  from  the  theory.  This  case 
will  be  studied  in  the  next  Section. 
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We  end  these  general  considerations  by  exhibiting  a  conservation 
law  intrinsic  to  the  general  formalism,  provided  the  system  is  autono¬ 
mous,  i.e.,  F  is  independent  of  t.  Given  any  autonomous  Lagrangian 
system,  the  Legendre  transform  of  the  Lagrangian  with  respect  to  the 
time  derivatives  of  the  dynamical  variables  (the  "energy")  is  time  inde¬ 
pendent.  A  straightforward  application  of  this  theorem  to  the  Lagrangi¬ 
an  appearing  in  eq.  (2.5),  viz.,  I,  =  l/2< (d^X-F) , (dtX-F)>  leads  to  the 

conservation  of  the  auasi-enerev 

Q  =  l/2<(dtX-F) , (etX+F)>.  (2.8) 

Unfortunately,  the  quasi-energy  is  not  a  positive  functional.  However, 
just  like  the  Lagrangian,  Q  vanishes  for  an  unperturbed  (6  =  0)  solu¬ 
tion  of  the  field  equations.  Therefore,  its  magnitude  characterizes  the 
randomness  of  the  response  of  the  system  under  the  influence  of  arbi¬ 
trarily  small  random  perturbations. 

3.  TURBULENCE  IN  AN  INCOMPRESSIBLE  VISCOUS  FLUID 

We  now  apply  the  general  formalism  outlined  in  the  previous  section 
to  the  turbulent  flow  of  a  viscous,  incompressible  fluid.  The  equations 
governing  the  flow  are: 

n*[v]  e  a  v*  +  (vsa  )v*  +  a*p  -  -  ef*,  (3.1) 

L  S 


where  v  stand  for  the  components  of  the  velocity,  p  is  the  pressure  and 
f  is  a  perturbing  Gaussian  random  force  with  unit  correlation  operator. 
In  order  to  apply  the  formalism  outlined  in  Sec.  2,  we  introduce  the 
momentum  canonically  conjugate  to  the  velocity  field,  w^ .  Then,  using 

eq.  (2.4),  the  action  entering  the  expression  of  the  characteristic 
functional  reads: 

A  =  /d3xdt{WjN* [v]  -  1/2  e2w*w^}.  (3.3) 

One  readily  convinces  oneself  that  the  part  of  the  action  linear  in 
is  invariant  under  the  following  scale  transformation: 

t  =  At';  X  =  A^2x'  ;  V  -  A  1/2V'  ;  p  =  A  *p'  ;  w  =  A  ^w',  (3.4) 

where  A  is  an  arbitrary  positive  scale  factor  and  v,  w  stand  for  the 
magnitudes  of  the  velocity  and  momentum  vectors.  (This  is  just  the 
usual  scaling  law  of  the  Navier-Stokes  equations;  in  particular,  the 
Reynolds  number  is  an  invariant  of  the  transformation.)  Therefore,  upon 
-4 

choosing  A  =  t  ,  the  coefficient  of  the  term  quadratic  in  the  canonical 
momentum  becomes  unity:  the  small  parameter  disappeared  from  the  theory 
altogether.  (More  precisely,  it  was  shifted  into  the  boundary  condi- 
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tions,  which,  however,  can  be  adjusted  accordingly.)  We  can  now  go 
directly  to  eq.  (2.5);  the  "Lagrangian"  form  of  the  action  reads: 


A  =  1/2  S d3xdt  NA[v]  N.  [v] . 


(3.5) 


This  action  has  to  be  supplemented  by  the  ghost  contribution,  r  , 
cf.,  eqs.  (2.6)  and  (2.7).  On  denoting  the  ghost  fields  by  £  and  q  , 
respectively,  we  have: 


r  =  /dt  d3x{£i(dt-yvt‘-vrc»r)7i 


*iV 


(3.6) 


In  order  to  develop  a  systematic  loop  expansion  to  the  functional  inte¬ 
gral  with  this  action,  we  need  to  determine  an  average  flow  which  is  a 

solution,  V  ,  of  the  unperturbed  Navier-Stokes  equations,  N  [V]  =  0  and 

of  the  subsidiary  condition,  d^V1  =  0-  The  solution  to  these  equations, 

together  with  £  =  q  =  0  (which  is  the  only  saddle  point  of  r),  determine 
a  saddle  point  of  the  integrand  in  the  functional  integral  expression  of 
the  generating  functional.  Let  us  suppose  that  such  a  solution  has  been 
found.  We  can  then  expand  the  action  around  the  saddle  point  noting 

that  the  fluctuations,  U*  =  v*  -  \T ,  have  to  be  divergence-free.  In 

order  to  satisfy  the  transversality  condition,  T  =  d^u1  =  0,  one  has  to 


insert  the  standard  Fadeev-Popov  factor,  6 (T) • Det(6T/6u) ,  under  the 
functional  integral.  The  determinant  is,  however,  independent  of  the 
integration  variables,  and  thus  no  ghost  has  to  be  introduced  for  this 
constraint.  It  is  further  a  straightforward  matter  to  show  that  the 
effect  of  the  delta  functional  is  merely  to  project  out  the  transverse 


part  of  the  propagator, 


by  means  of  the  projector, 


T.  . 
ij 


-2^  ^ 

-  V  d  ,d  .  . 
1  J 


i  i  i  i 

After  inserting  v  =  V  +  u  into  the  action  and  using  N  [V]  =  0, 

the  action  splits  up  as  follows: 


A  =  A[V]  -t-  AQ(u]  +  A1tu],  r  =  rQ[£  ,rj  ]  +  rj[£  ,f7,u]  ,  (3.7) 

where  the  subscripts  0  and  1  stand  for  the  quadratic  ("kinetic")  and 
higher  order  ("interaction")  parts  of  the  action  describing  the  fluctua¬ 
tions,  respectively.  The  general  expression  of  the  various  pieces  can 
be  best  described  as  follows.  Introduce  the  following  quantities: 


1.  =  (3  -  wv2  +  V1^  )u.  +  V.  ur, 

i  t  r  i  i  f  r 


Lik  -  <st  -  *  vrV6ik +  vi,k 


We  then  have 


(3.8) 
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Aq[u]  =  1/2/dt  d3x[l.l1  + /i"2^1)2]. 

Aj[u]  =  1/2/dt  d3x{2(l1uku1  R)  +  (u*uk  1)(uruk,r)}, 
rQ  =  /dt  d3xt1L1Rf7k, 

ri  =  /dt  d3xt  ^ (urdr6^k  +  Uj  y)qR  .  (3.9) 

From  here  one  Immediately  reads  off  the  general  Feynman  rules  needed  for 
the  computation  of  the  non-Gaussian  part  of  the  distribution.  In  par¬ 
ticular,  the  ghost  propagator  is  given  by  L  1  whereas  the  propagator  of 

the  velocity  fluctuations  is  given  by  (L*L)  *.  Clearly,  in  order  to 
proceed  any  further,  the  average  flow  has  to  be  specified.  In  the  next 
section,  this  procedure  is  illustrated  for  the  important  example  of  a 
flow  which  is,  on  the  average,  Poiseuille. 

4.  FLUCTUATIONS  IN  A  FLOW  BETWEEN  PARALLEL  PLANES 

We  consider  a  flow  taking  place  between  two  parallel  planes;  the 
pressure  gradient  is  stationary  and  homogeneous;  its  direction  is  paral¬ 
lel  to  the  planes.  We  assume  that  the  average  flow  is  given  by 
Poiseuille's  solution  of  the  Navier-Stokes  equations,  viz.,  in  a 
suitably  chosen  Cartesian  system  of  coordinates,  the  nonvanishing 
component  of  the  average  velocity  is 

V  =  1/2  U ( 1  -  x2/a2),  (4.1) 

z 

where  U  is  related  to  the  pressure  gradient  and  kinematic  viscosity  in  a 
well-known  way,  and  "a"  stands  for  the  half-distance  between  the  planes 
confining  the  flow.  The  coordinate  axes  x  and  z  are  perpendicular  and 
parallel  to  the  planes,  respectively,  whereas  the  y-axis  is  perpendic¬ 
ular  to  the  (x,z)  plane.  Our  objective  is  to  work  out  the  expression 
for  the  action  with  this  average  flow,  and  then  carry  out  the  inversion 

of  L  and  L*L  in  order  to  obtain  the  Feynman  rules  which  ultimately  de¬ 
termine  the  correlation  functions.  Even  with  this  simple  average  flow, 
this  is  a  rather  tedious  task  and  the  result  cannot  be  obtained  in  a 
closed  form.  In  order  to  proceed,  we  rescale  all  quantities  in  order  to 
measure  them  in  their  natural  units.  To  this  end,  the  unit  of  length  is 
chosen  to  be  a,  the  unit  of  time  is  the  characteristic  advection  time, 

a/U,  and  velocities  are  measured  in  units  of  U.  In  addition,  the  fluid 

* 

is  characterized  by  an  intrinsic  diffusion  velocity,  u  =  iVa.  The 

* 

Reynolds  number  is  defined  to  be  R  =  U/u  .  We  can  choose  the  units  in 
which  the  correlation  operator  of  the  stirring  is  measured  such  that  it 

*  3  g 

is  proportional  to  (u  )  a  .  In  this  way,  if  wo  go  over  to  dimensionless 
quantities,  the  quadratic  part  of  the  action,  A^,  will  be  proportional 
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3 

to  R  .  In  these  units  the  average  flow  is  obviously  given  by  V  = 

z 

2  k 

l/2(l-x  ).  Next,  define,  as  in  the  preceding  section,  lj  =  L^u  ,  where 

k 

u  is  the  velocity  fluctuation  in  its  natural  units  and  is  the 
operator  (again  using  natural  units): 

Lik  “  (\  -  R'v  *  vrV5ik  *  vl'k-  <4-2* 

In  this  way  we  get, 

Aq  =  R3/dtd3xl1l . .  (4.3) 

It  is  now  obvious  that  at  large  Reynolds  numbers  the  action  is  propor- 

3 

tional  to  the  large  number,  R  .  In  the  same  way  as  in  quantum  field 
theory  the  semiclassical  limit  (h  -*  0)  is  obtained  by  rescaling  the 
fields  with  an  appropriate  power  of  Planck's  constant;  here  we  remove 

i  _ 3 /2  x 

the  large  number  in  front  of  AQ  by  rescaling,  u  ->  R  ,  u  .  In  this 

way,  the  quadratic  part  of  the  action  becomes  0(1)  in  R,  whereas  the 
non-Gaussian  part  carries  negative;  powers  of  R;  thus  it  is  relatively 
unimportant  for  large  Reynolds  numbers.  (In  fact,  one  can  define  a 
systematic  expansion  in  inverse  powers  of  R;  here  we  restrict  ourselves 
to  the  leading  term  of  that  expansion.)  It  is  also  easily  seen  that  in 
the  leading  approximation  to  the  distribution,  the  ghost  fields  decouple 
from  the  velocity  fluctuations.  Thus,  in  the  leading  approximation,  one 
has  a  Gaussian  distribution  of  fluctuations  for  R  >>  1.  In  order  to 
compute  the  correlation  function  of  the  fluctuations,  we  now  have  to 

invert  the  operator  i/l.  This  is  best  done  in  a  mixed  representation: 

we  Fourier  transform  in  the  variables  t,  y,  z,  denoting  the  conjugate 

variables  by  w,  k  k  ,  respectively.  One  is  tempted  to  omit  the  term  of 
y  z 

0(R  )  in  the  expression  of  L  altogether.  However,  in  this  way,  the 

viscous  effects  would  be  totally  neglected.  Instead,  at  least  for  the 

- 1  2  2 

purpose  of  illustration,  we  omit  the  term  R  d  /dx  from  eq.  (4.2),  thus 

6 

making  the  distribution  independent-valued  in  x.  This  is  expected  to 
be  a  reasonably  good  zeroth  approximation  outside  the  boundary  layer 

-1/2 

which,  in  natural  units,  is  approximately  of  thickness  R  .  In  this 
approximation,  the  computation  of  the  velocity  correlation  function  (the 
propagator  of  the  fluctuation)  is  a  straightforward,  although  somewhat 
tedious,  task:  the  technical  complication  arises  from  the  application 

of  the  nonlocal  projector,  T^.,  on  (L*L)  * .  If  we  are  to  remain  in  the 

leading  approximation,  however,  we  do  not  need  the  full  propagator.  The 
correlation  function  of  vorticities  is  easier  to  compute:  one  just  has 


6 

I.  M.  Gel'fand  and  N.Ya.  Vilenkin,  Generalized  Functions.  Vo] .  4, 
Academic  Press,  New  York,  1964,  Chapter  III. 
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to  take  the  curl  of  a  matrix  element  of  the  operator  (L  L)  in  both  of 
Its  variables.  (Obviously,  for  an  incompressible  fluid,  the  vorticity 
correlation  function  contains  the  same  physical  information  as  the  full 
correlation  function.)  Denoting  the  components  of  the  vorticity  cor¬ 
relation  tensor  by  the  result  of  a  straightforward  computation  can 

be  described  as  follows. 


i)  From  the  reality  and  symmetry  of  the  vorticity  correlation 
tensor  in  coordinate  space,  one  readily  derives  a  Hermiticity  relation¬ 
ship  in  the  mixed  representation  used  here,  viz., 

* 

nij(~ky’_1V~W:X,X ')  =  nji(ky’kz’w  :x'  ,x>  *  (4.4) 

ii)  Taking  this  into  account,  we  list  a  set  of  independent  vorti¬ 
city  correlation  functions  as  follows. 


Let  us  define  the  quantity,  B  =  (i[w+l/2(l-x2)k  ]  +  R  1(k2  +  k2)}. 

z  y  z 

Then, 

Q  =  — I B I  2 [k2 (1  +  X  )  +  k2 ]5 (x  -  x' ) , 
xx  11  1  y '  | B ) 2  Z 


2  2 
~{|b|~2(1  +  — — (x  -  X')} 
dx  I  B( 


-  ikz  ^{XR~V2  +  k2 )  |  R  j  46(x  -  x')} 


5  ( x  -  x '  )  , 


n 


zz 


2 

— -g* [  | B j  26(x  -  x’)]  -  ky|B|  25(x  -  x  1  )  , 
dx  y 


a 

yx 


■  -  V1  §3r‘|Bf*<i 


+  -^-^)G(x  -  x’ )  ]  +  xk  6{x  -  x 1 )  > , 

|B|2  2  |B|4 


n 

zx 


lkz  l?[lBr26(x  -  X')] 


xk 


2  B 


6(x 


(l  =  k  (k  B  26(x  -  x')  +  i  |— [  — 

zy  y  z  a  x L - B -  4 


x6(x  -  x ' ) ] } 


(4.5) 


Not  surprisingly,  the  correlation  functions  are  proportional  to  6-func¬ 
tions  and  their  derivatives  in  x:  this  is  to  be  expected  for  indepen¬ 
dent-valued  distributions,  cf . ,  Gel'fand  and  Vilenkin,  loc.  cit. 
Accordingly,  one  has  to  fold  the  correlation  functions  quoted  with  the 
efficiency  functions  of  the  detectors  in  order  to  get  physically  mean¬ 
ingful  expressions. 
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5.  DISCUSSION 


The  formalism  developed  in  this  paper  allows  one  to  develop  a 
systematic  approximation  procedure  to  the  theory  of  stochastic  systems 
originally  developed  in  refs.  1  and  2.  In  particular,  as  it  was  demon¬ 
strated  in  the  preceding  section,  one  can  obtain  useful  analytical 
results  in  the  form  of  an  asymptotic  expansion  in  inverse  powers  of  R 
in,  at  least  some,  turbulent  flows.  However,  the  usefulness  of  this 
technique  goes  beyond  the  large  R  expansion.  In  its  present  form,  the 
formalism  appears  to  be  much  more  amenable  to  the  application  of 
standard,  non-perturbative  techniques  extensively  used  in  quantum  field 
theory,  such  as  a  renormalization  group  analysis  of  the  various  correla¬ 
tion  functions.  The  technical  reason  for  this  is  that  the  present, 
Lagrangian  formulation  of  the  functional  integral  is  better  adapted  to 
the  application  of  such  techniques  than  the  essentially  canonical 
formulation  in  refs.  1  and  2.  In  those  papers,  one  was  forced  to  use 
retarded,  rather  than  Feynman  propagators:  this  circumstance  alone 
rendered  the  formalism  so  complicated  that  an  efficient  implementation 
of  modern  techniques  was  virtually  impossible. 

In  addition,  it  appears  that  using  the  proposed  formalism  one  can 
gain  some  interesting  insights  into  the  theoretical  description  of 
chaotic  behavior  in  general  and  the  theory  of  turbulence  in  particular. 

7 

It  has  been  emphasized  in  the  literature  that  chaotic  behavior  has  many 
features  which  are  similar  to  phase  transitions  or  dynamical  symmetry 
breaking.  In  this  respect,  it  is  particularly  pleasing  to  notice  that  a 
consistent  statistical  description  of  chaotic  phenomena  appears  to 
require  the  presence  of  an  arbitrarily  weak  but  nonvanishing  "disorder¬ 
ing  field,"  just  as  a  consistent  theory  of  dynamical  symmetry  breaking 

g 

requires  the  presence  of  a  weak  "ordering  field." 


7 

Hao  Bai-Lin,  Chaos .  World  Scientific,  Singapore,  1984,  Chapter  1. 

8 

G.  Domokos  and  P.  Suranyi,  "Spontaneous  Symmetry  Breaking  in  Quantum 
Field  Theory,"  Sov.  Jour.  Nuc.  Phvs . ,  Vol .  2,  pp.  361-367 ,  1966. 
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